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SUMMARY 


The model for prediction of low-frequency stability limits for bipropellant rocket 
engines is modified to include a discrete vaporization time for each propellant, plus a 
mixing and reaction time common to both propellants. With this modified model as a 
hypothesis, stability boundaries are calculated for an engine burning liquid oxygen and 
gaseous hydrogen. The results obtained are, for certain combinations of vaporization 
and mixing times, converse to those which would be obtained by a lumped -dead -time 
model such as has been used. For these cases in particular, stability can be improved 
by decreasing the hydrogen-injector pressure drop, or by increasing the oxygen- 
vaporization time. Until now, these two concepts have been diametrically opposed to 
design criteria for stability improvement. The limited data that exist qualitatively sup- 
port the modified model. 

For simplicity, factors such as a complex feed system and sensitivity of the dead 
times to chamber pressure and injection velocity are not included in this analysis^ J 


INTRODUCTION 


Low-frequency stability of liquid -propellant rocket engines has been the subject of 
many analyses during the past two decades. The analytical models currently in use to 
predict stability limits have evolved from the contributions of many researchers. Von 
Karman (ref. 1) introduced the concept of a combustion delay (dead time) as an explana- 
tion of chamber pressure oscillations, Gunder and Friant (ref. 2) added feed system 
inertance, and Summerfield (ref. 1) introduced combustion chamber dynamics to the 
model. Crocco and Cheng (ref. 3) refined the model by letting the dead time vary with 
chamber pressure, and Hurrell (ref. 4) added injection velocity effects. 

In previous analyses, the vaporization time of the propellant having the longer drop 


lifetime was calculated in the manner of Priem and Heidmann (ref. 5). This vaporiza- 
tion time, generally referred to as dead time, was then applied to both propellants. In 
rocket engines where the propellants have radically different drop lifetimes or where one 
of the propellants is introduced as a gas, stability models with a common dead time 
yield questionable results. 

This report advances a stability model in which each propellant is acted on by a dis- 
crete dead time. This concept is used in the development of the equations describing the 
stability boundaries, and typical analytical results are presented. Since it is the pur- 
pose of this report to stress the importance of discrete dead times, neither the effects of 
the feed system nor the pressure and velocity sensitivity of the dead times will be con- 
sidered. 


ANALYSIS 

Combustion Process 

The combustion process in a liquid-propellant rocket engine is still not clearly 
understood, although what takes place is generally known. Initially, the propellants are 
injected into the combustion chamber; their flow rates are determined by the conditions 
upstream and downstream of the injector element and the element geometry. The pro- 
pellants are atomized, vaporized, and mixed, and then they react to produce hot gases. 
These processes are gradual and continuous. For mathematical expendiency, however, 
they are treated in a discontinuous manner; that is, the gradual evolution of several 
small elements of propellants into a particle of burned gas is replaced by a sudden con- 
version. The time interval between injection and sudden conversion is called dead time. 
The manner in which dead time is applied to propellants will be critically reviewed in 
this report. 

As noted in reference 1, the time rate of change of chamber pressure is determined 
by the rate of generation of burned gases and by the rate of depletion of the burned gases 
through the nozzle. The characteristic time associated with this process is the well- 
known gas -residence time 6 . 


New Model 

In previous analyses, the dead times for all the processes were lumped together and 
applied equally to both propellants. While it is true that mixing and reaction times are 
common, vaporization times are associated with individual propellants and should be 
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(a) Single-dead-time model. 
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(b) Postulated model. 

Figure 1. - Block diagrams for stability limit models for bipropellant rocket engines. 


treated accordingly. These concepts are illustrated as block diagrams in figure 1. For 
comparison, the formerly used single -dead -time model is presented in figure 1(a), and 
the postulated model with dead times acting on their respective flows is presented in 
figure 1(b). In this figure, the feed system is assumed to be completely decoupled. It is 
important to note the discreteness of the dead times, that is, the assigning of different 
dead times to each flow loop. This is the primary departure from the models used in the 
past. The following equations describe this model. (Symbols are defined in the ap- 
pendix. ) 

By definition 


P x ,o(t)-P C < t > = 4P I,o<‘) 

(1) 

P T f(t> - P c (t) = AP Ijf (t) 

(2) 


The flow through the injectors can be expressed by 
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v t,=K °h , ,o ( ‘> 


L 

where the K’s and the exponents are assumed to be constant, 

Wf, v ^ = W f>i (t-a v>f ) 
W,,(t)=W 0 >v (t-a m ) 


o,b 




(3) 

(4) 


(5) 

( 6 ) 

(7) 

( 8 ) 


By conservation of mass, 


ih<*>] = * 0 ,b<‘) * , v< t) - 

If the burned products behave as a perfect gas, 


(9) 


RT C 

F cW - -j- w s (t) 


( 10 ) 


For a choked nozzle 


w N (t) = 


P c (t)Ag 

C*(t) 


(ID 


After linearization and/ or Laplace transformation, these equations become 


dW o i 

6 ^o,i (s) = dAPj D [ 6 P T,o ( s) ‘ 6P C (s) ] 


dW f.i 


6W l,i (S,= 35r[ 6P T,£ (S) - 6P C (s) ] 
1> ® 


( 12 ) 

(13) 
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(14) 


6W 0jV (s) = 6W 0ji (s)e'V oS 


(15) 

(16) 

(17) 

(18) 

The gas residence time 6 is defined as 

o 

c Vp P*T * 

0 = — = -^=- (19) 

g RT c Ag RT c g 


6W. 


-a ,s 


W fv (s) = 8W f;i (s)e v ’ f 


6W 0ib (s) - 8W 0>v ( s )e 


-o s 
m 


8P c (s) = 


5W f,b (s) = 6Wf > v 

'dP, 


(s)e 


-cr s 
m 


9gS+ i 


L' 3W o- 




and is assumed to be constant. 

Equation (10) to (17) can be combined to yield 



( 21 ) 
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The denominator of equation (19) is, when set equal to zero, the characteristic 
equation 


- ct s 
m 


0=1 + 


0„s + 1 
g 




k 0W Ad AP 


I,o/ 


^ W f/\ dAP I,i 


( 22 ) 


(see, e.g., ref. 6). Evaluating the derivatives results in 


dW o a-1 W Q 

= alC AP T = a 


d AP- 


I,o 


V LIT -r — CL 

o 1, 0 AP 


1,0 


(23) 


Since K Q = (W^/AP^), and 


dW f W f 

= b 


dAP I,f ^I,i 


(24) 


Ppi = “7 (W i + W f . ) 
C Ag °,b f,b 7 


(25) 


3P 


C 1 


9W o,b 


C* + (w„ h + w. J 


ac” 


°> b f,b ; gw 


o,b 


(26) 


ac* ac* dR ac* 1 


aw n . 3R dW^ h dR W f h 

o,b o,b f,b 


(27) 


where R = W Q and 


3P, 


3W 0)b Ag 


C* + (R + 1) 


ac* 

3R 


(28) 


In like manner, 


3P 


C 1 


3W f Ag 


c* - R(R + l) 


ac* 

3R 


( 29 ) 


Substituting equations (23), (24), (28), and (29) into equation (22) gives 
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-CT_S f „tt7 


m 


aW, 


Ag 0 g s + nAP I o 


C* + (R + 1) — 

-a s bW. 
e v -° ♦ f 

C* - R(R + 1) 

e ~ CT v, f s J 

3R_ 
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9R_ 
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(30) 


Equation (30) rewritten in terms of APj/Pq is 



Equation (31) may be reduced to a single -dead-time case by setting A y Q equal to 

and combining with A so that ’ 

m 




(32) 


RESULTS 

Equation (31) may be implemented on a digital computer to solve for engine stability. 
For example, consider an engine running on gaseous hydrogen and liquid oxygen. Since 
the hydrogen enters the chamber as a gas, fuel vaporization time is zero. In figure 2, 
stability boundaries are presented for an oxidant vaporization time of 2 milliseconds and 
mixing times varying between 0 and 4 milliseconds. A point on the stability boundary is 
determined by those values of APj j/Pq and APj Q /P^, which yield one or more conju- 
gate pairs of roots of equation (29) which lie on the imaginary axis (s = a + ja>; a - 0), 
with the condition that no roots exist to the right of the imaginary axis (s = a + jcu; 
a > 0). The stability boundaries are presented as loci of these points. 

In figure 2, the stable operating region lies above and/or to the right of the boundary 
in each case. It is seen that, as mixing time is decreased, the area of stable operation 
increases. Furthermore, the shape of the curves becomes such that a transition from 
unstable to stable operation can be achieved by decreasing the fuel pressure drop. 
Although contrary to the usual practice of increasing pressure drops to stabilize an 
engine, this concept is in agreement with experimental experience on an engine of this 
type (ref. 7) where this exact behavior was observed. 
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Figure 2. - Effect of mixing time on stability boundary. Vaporization time, 2.0 milliseconds. 

The frequency at which the oscillations occur changes along the curve, although this 
fact is not indicated in figure 2. The discontinuities in the curves (for mixing times be- 
tween 0. 1 and 0. 8 msec) correspond to a flip in frequency from one mode to another. At 
the inflection points, two pairs of roots exist on the imaginary axis and oscillations occur 
at both the associated frequencies. 

Study of the behavior predicted by this model indicates that stability may be 
achieved by manipulation of the dead times. This point is illustrated in figure 3. Mixing 
dead time is held constant at 1. 0 millisecond, and stability boundaries are plotted for 
vaporization dead times ranging from 0. 5 to 2.0 milliseconds. Again, the stable region 
lies above and/or to the right of each curve. If the operating point indicated by the 
symbol were chosen, a transition from unstable to stable operation could be realized by 
increasing the oxidant vaporization time from some lower value to 1. 5 milliseconds. 

This behavior is contrary to that which would be predicted by a single -dead -time model, 
where any increase in dead time adversely affects stability. It should be pointed out, 
however, that further increases in dead time are destabilizing. 

The behavior predicted by the single -dead -time model, which is described by equa- 
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Figure 3. - Effect of vaporization time on 
stability boundary. Mixing time, 1.0 
millisecond. 
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Figure 4. - Stability boundaries calcu- 
lated with single-dead-time model. 


tion (32), is illustrated in figure 4 for mixing times ranging from 1.0 to 3. 0 milli- 
seconds. For this model, any increase in dead time, or decrease in pressure drop, is 
destabilizing. 


CONCLUDING REMARKS 

This report has advanced the concept of separate dead times and demonstrated the 
effects of using this concept in chugging analyses. Since the presentation of results 
would have been encumbered by variables such as a complex feed system and the sensi- 
tivity of the dead times to pressure and injection velocity, these effects have not been 
included. 

There are limited experimental data to qualitatively substantiate the behavior pre- 
dicated by this model. The model does, however, provide an explanation for the data in 
reference 7, where chugging was eliminated by decreasing the fuel injector pressure drop. 
This result is inexplicable when a single-dead-time model is used. 

Reference 8 indicates that, for hydrogen-oxygen engines, a low oxygen -injection 
velocity is desirable from a screeching standpoint. This implies a long oxygen- 
vaporization time, which has previously been deemed undesirable from a chugging point 
of view. The model advanced herein indicates that a long vaporization time may not be 
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undesirable; hence, an injector with a low propensity to screech may not compromise 
low-frequency stability. 

If the performance predicted by this model is corroborated as additional data are ob- 
tained, currently used design criteria such as maintaining a minimum AP/P^ and 
striving for low vaporization times should be critically reviewed. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 5, 1965. 
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APPENDIX 


SYMBOLS 


A throat area, in. 

a exponent of oxidant injector 

b exponent of fuel injector 

C* characteristic exhaust velocity, 

in. /sec 

g gravity, in./sec^ 

j V=i 

K injector resistance constant 

L* characteristic chamber length, in. 

o 

P pressure, lb/in. abs 

R mixture ratio, 

(ft gas constant, in./°R 

s Laplace operator, 1/sec 

T temperature, °R 

t time, sec 

3 

V volume, in. 

W weight, lb 

W propellant flow rate, lb/ sec 


a real part of Laplace operator, s 

5 perturbation 

gas residence time, sec 

0 dead time, sec 

co imaginary part of Laplace 

operator, s 

Subscripts: 

b burned 

C chamber 

f fuel 

1 injector 

i injected 

m mixing, reaction, and heating 

N nozzle 

o oxidant 

s stored 

T tank 

v vaporized, vaporization 
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